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FOREWORD 


This  report  describes  research  performed  by  Universal  Energy 
Systems,  Inc.,  Dayton,  Ohio.  The  work  was  conducted  under  Contract 
F33615-77-C-31 13 ,  "Investigations  of  Methods  of  Plasma  Excitation," 

Task  III,  "Advanced  Optical  Diagnostics  of  Plasmas." 

The  work  reported  herein  was  performed  at  in-house  facilities  of 
the  Air  Force  Aero  Propulsion  Laboratory,  Air  Force  Wright  Aeronautical 
Laboratories,  Wright-Patterson  Air  Force  Base,  Ohio.  The  experimental 
work  and  analysis  described  was  performed  by  C.  A.  DeJoseph.  As  Mr. 
DeJoseph  was  not  employed  by  Universal  Energy  Systems  at  the  time  of  the 
completion  of  the  contract,  this  section  of  the  final  report  was  prepared 
by  V.  E.  Merchant  by  combining  and  editing  quarterly  reports  previously 
written  by  Mr.  DeJoseph. 
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INTRODUCTION 


Advanced  optical  techniques  have  been  applied  to  several  problems  in 
plasmas  diagnostics.  A  technique  has  been  developed  for  determining  the 
rotational  temperature  and  vibrational  population  distribution  of  non- 
homonuclear  molecules  under  non-equilibrium  conditions.  The  technique 
was  applied  to  carbon  monoxide  under  the  conditions  present  in  an 
electrical  discharge.  The  second  problem  area  was  identifying  infrared 
absorption  spectra  of  impurities  that  might  be  present  in  C02  laser  dischar¬ 
ges.  Spectral  regions  free  from  interference  from  the  primary  constituents 
and  common  impurities  such  as  water  vapor  were  identified  for  several  of  the 
nitrogen  oxides,  and  the  detection  sensitivities  were  determined.  Finally 
this  information  was  used  in  a  study  of  impurity  buildup  during  discharges 
in  an  electron-beam  excited  closed  cycle  (EBCC)  system  and  in  a  TEA  laser. 

Succeeding  sections  of  this  report  give  more  information  about  the 
developments  in  these  main  areas  of  concentration.  In  addition  several 
minor  problems  were  addressed.  The  infrared  transmittance  of  Freon  1301 
in  the  10.6  pm  spectral  region  was  investigated.  Infrared  emission  spectra 
were  obtained  from  He-Ne  and  He-Xe  hollow  cathode  discharges.  No  further 
information  about  these  minor  problem  areas  will  be  given  here. 

Throughout  the  duration  of  the  contract,  several  improvements  were 
made  to  the  Fourier  Transform  Spectrometer  (FTS).  A  variable  path  corrosion 
resistant  absorption  cell  and  associated  gas  handling  hardware  were  added 
to  the  system  for  use  with  corrosive  gases.  The  system  constructed  for  the 
trace  species  absorption  measurements  is  centered  around  a  long-path  (20  meters) 
White-type  infrared  absorption  cell  with  optics  for  use  with  the  FTS  system. 


In  addition,  a  stainless  steel,  high  vacuum  gas  handling  system  is  used  for 
evacuating  the  cell  and  preparing  the  synthetic  mixes.  Figure  1  shows 
the  layout  of  the  optics  system  for  the  White  cell.  The  glowbar  source  is 
imaged  1:1  onto  the  adjustable  aperture  by  mirror  #1  and  then  collimated  by 
mirror  P2.  Mirror  P3  focuses  the  radiation  onto  the  entrance  window  of  the 
White  cell  and  the  radiation  leaving  the  cell  is  re-collimated  by  mirror  H. 
Radiation  leaving  mirror  PA  then  passes  to  the  FTS  System  where  the  spectrum 
is  measured.  Figure  2  shows  the  gas  handling  system  for  the  White  cell. 

High  vacuum  evacuation  of  the  system  is  by  the  liquid  nitrogen  trapped  diffusion 
pump.  Base  pressure  for  the  system  is  the  mid  10" 7  torr  range.  The  system 
can  be  used  in  two  ways.  Samples  of  gas  from  the  EBCC  loop  can  be  entered 
into  the  system  via  the  sample  port  or  absorption  spectra  of  various  gas  mixes 
can  be  obtained  by  entering  these  into  the  system  through  the  gas  manifold. 

Gas  pressure  was  measured  with  the  two  MKS  Baratrons  which  gave  an  increase 
in  dynamic  range  of  pressure  measurement  over  a  single  unit.  Gas  mixes 
were  prepared  under  static  pressure  conditions  and  likewise  the  absorption 
measurements  were  made  under  static  conditions.  After  a  gas  fill,  the 
White  cell  was  "roughed  down"  using  forepump  P 2  which  is  a  high  throughput 
Leybold  Heraeus  model  DK-45.  Use  of  this  pump  greatly  facilitated  pump  down 
of  the  cell  from  atmospheric  pressure  while  allowing  forepump  #1  to  remain 
small  for  reduced  backstreaming  to  the  diffusion  pump. 

The  spectral  range  of  study  was  originally  limited  to  the  range  of 
5.9  urn  to  1.7  pm  (1700-6000  cm-1)-  This  is  the  range  over  which  the  original 
PbSe  detector  used  in  the  FTS  System  was  sensitive.  A  new  "sandwich" 
detector  was  obtained  which  allowed  the  range  to  be  extended  from  14  ym 
to  1.7  ym  (700-6000  cm-1).  This  detector  is  composed  of  a  LN2  cooled 
HgCdTe  detector  "sandwiched"  behind  a  LN2  cooled  InSb  detector.  The  detector 
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FIGURE  1  WHITE-CELL  OPTICS 


can  he  wired  so  that  the  outputs  are  added  or  come  out  separately,  with 
each  configuration  having  its  own  advantages  and  disadvantages.  For  maximum 
signal-to-noise,  the  outputs  are  taken  separately  and  either  part  of  the 
composite  detector  is  selected  by  the  FTS  System  controller.  Thus,  in  its 
present  configuration  the  FTS  System  can  gather  data  in  the  range  of  the 
HgCdTe  detector,  ie.  14  pm  to  5.4  pm  (700-1850  cm'1)  or  the  InSb  detector, 
ie.  5.7  pm  to  1.7  pm  (1750-6000  cm-1).  Results  shown  in  this  report  were 
primarily  obtained  with  the  PbSe  detector.  Most  measurements  were  taken  in 
the  presence  of  laboratory  air,  so  it  was  necessary  to  contend  with  parasitic 
absorption  by  water  vapor  and  carbon  dioxide  in  the  air.  Figure  3  shows 
a  low  resolution  spectrum  of  the  glowbar  source  as  seen  through  the  White 
cell  in  the  InSb  (or  roughly  the  PbSe)  region  of  the  spectrum  in  lab  air. 

The  regions  showing  absorption  due  to  H-0  and  CO-  must  be  taken  into  account 
when  deciding  where  to  look  for  absorption  by  the  minor  species.  Figure 
3  will  serve  to  locate  absorption  features  shown  in  more  detailed  spectra 
later  in  the  report. 

A  significant  change  to  the  absorption  cell-FTS  optics  was  the  installing 
of  a  germanium  on  potassium  bromide  (KBr)  beamspl itter  in  the  FTS  system.  With 
this  beamsplitter  and  the  "sandwich"  detector  spectral  coverage  of  the  FTS 
system  has  been  extended  to  the  range  14  pm  to  2  pm.  Because  of  the  way  the 
"sandwich"  detector  is  wired,  the  spectral  range  must  be  covered  in  two 
separate  data  acquisitions.  Data  is  acquired  over  the  spectral  range  14  pm 
to  5.4  pm  (700-1850  cm-1)  or  5.7  pm  to  2  pm  (1750-5000  cm"1).  Thus,  analysis 
of  a  gas  sample  requires  two  data  acquisitions  each  requiring  about  twenty 
minutes  to  complete.  Replacing  the  iron  oxide  on  calcium  flouride  (CaF:) 
beamsplitter  with  the  KBr  extends  the  spectral  region  from  about  1100  cm'1 

to  700  cm-1  on  the  low  frequency  side  while  shortening  the  region  from 
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6000  cnr 1  to  5000  cm-1  on  the  high  frequency  side.  However,  due  to  the 
importance  of  the  10  pm  spectral  region  in  some  experiments,  the  increased 
coverage  in  the  low  frequency  part  of  the  spectrum  is  a  more  desirable 
configuration. 

An  investigation  was  made  of  the  factors  limiting  the  precision  of 
wavelength  determination  with  the  Fourier  Transform  Spectrometer.  The 
precision  with  which  a  line  position  can  be  determined  is  dependent  on  a 
number  of  experimental  parameters.  One  parameter  that  is  fundamental  to  the 
best  achievable  precision  is  the  resolution  of  the  instrument  12] .  If  the 
only  factor  affecting  the  precision  of  the  measurement  is  the  experimenters 
ability  to  locate  the  center  of  the  line,  then  the  best  precision  achievable 
is  likely  to  be  between  one-twentieth  and  one-thirtieth  of  the  resolution  [2], 
The  AFAPL  FTS  system  possesses  a  maximum  resolution  of  0.6  cm-1  (apodized) 
and  a  theoretical  resolution  of  .037  cm-1  (unapodized).  Therefore,  the 
maximum  precision  with  which  a  line  position  can  be  determined  with  this 
instrument  should  be  on  the  order  of  ±.001  and  ±.003  cm'1.  A  study  was 
undertaken  to  determine  how  closely  this  maximum  could  be  approached.  Wave- 
number  precision  of  the  recovered  spectrum  from  the  FTS  system  is  directly 
dependent  on  the  precision  with  which  the  moving  mirror  position  is  measured 
in  the  interferogram.  The  mirror  position  is  determined  by  monitoring  the 
interference  fringes  produced  by  a  helium-neon  (HeNe)  laser,  so  the  wave- 
number  scale  in  the  recovered  spectrum  is  reasonably  precise  due  to  this 
technique  of  mirror  position  measurement.  This  "automatic"  calibration 
of  the  FTS-produced  spectrum  is  sometimes  referred  to  as  Conne's  advantage. 

A  number  of  factors  effect  the  best  achievable  precision  of  this 
"automatic"  calibration  scheme  so  an  attempt  was  made  to  identify  and  eliminate 
the  largest  of  these.  One  factor  effecting  the  precision  of  the  calculated 


wavenumbers  which  is  well  known  is  the  degree  to  which  radiation  from  the 
source  is  collimated  [2].  The  reason  for  this  is  that  rays  entering  the 
interferometer  off-axis  travel  a  greater  distance  than  those  entering  on-axis, 
thus  the  computed  spectrum  for  the  off-axis  rays  differs  from  the  spectrum 
computed  for  on-axis  rays.  As  Bell  has  shown  [2],  if  the  source  subtends 
a  solid  angle  n,  the  computed  wavenumbers  a are  related  to  the  true  wave- 
numbers  o  by; 

°c  =  °o  ^  (1) 

Thus,  a  measured  line  position  will  be  lower  in  wavenumber  than  the  true 
line  position.  Using  the  small  angle  approximation  for  equation  (1) 
can  be  written, 

°c  =  °0  *  °‘2/4) 

where  a  is  the  angle  of  deviation  of  rays  entering  the  interferometer. 

For  example,  if  one  requires  the  wavenumbers  to  be  within  +.01  cm-1  at 
4000  cm"1,  m  must  be  3.2  mrad.  or  11  minutes  of  arc.  Thus  it  is  necessary 
when  doing  precision  line  position  measurements  to  precisely  control  the 
degree  of  collimation  of  the  source. 

A  second  source  of  error  results  from  the  apparatus  operating  in  air 
and  the  dispersion  in  the  index  of  refraction  of  the  air.  Since  the  mirror 
position  is  measured  with  the  Hefle  laser,  the  wavelength  of  the  laser  serves 
the  same  purpose  as  a  known  line  in  a  conventional  spectrometer.  Because 
the  known  line  in  the  FTS  system  is  quite  removed  in  wavenumber  from  the  un¬ 
known  lines,  the  vacuum  correction  for  the  FleNe  laser  is  not  the  proper 
correction  for  the  IR  region  of  the  spectrum.  To  see  this  source  of  error 
more  clearly,  consider  the  fundamental  equations  of  FT-IR.  Let  the  fluctuating 
part  of  the  interferogram  as  a  function  of  path  difference  x  be  I (x ) ,  and 
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the  spectral  distribution  function  be  B(o).  Then  for  an  ideal  interferometer, 
I (x)  and  B(o)  are  related  by,  , 

I(x)  =  2[  B(o )COS(2na  x)d a  (2) 


0 

,oo 


B(a)  =  2 


I(x)C0S(2no  x)dx 


(3) 


where  o  is  in  cm-1,  and  the  constant  2  is  somewhat  arbitrary.  Assume  that 
a  source  possesses  a  spectral  distribution  function  (in  air)  of  B  (a  ).  The 
measured  interferogram  from  the  source,  I  (X  ),  is  then  given  by. 


Im(Xj  =  2 
m  m 


B  (o  )  COS  2tto  X  do  . 
o  o  o  m  o 


(4) 


The  calculated  spectrum  Bc(oc)  is  given  by, 


Bc(oc)  ■  2 


I  (X  )  COS  2 71-  X  dX  . 
mm  c  c  c 


(5) 


Substituting  equation  (4)  into  (5)  gives, 

Bc(ac>  =  41"[  f  W  C0S  2’“oXmdool  C0S  2™cXcdXc-  <6> 

0  Jo  -1 

Here  the  subscripts  'c*  and  'm'  refer  to  quantities  either  used  in  or  the 
results  of  calculations,  and  quantities  actually  measured.  The  measured 
interferogram,  I  ,  is  acquired  by  monitoring  the  fringes  of  the  HeNe  laser 
and  sampling  at  every  j1  zero-crossing  of  the  laser  interferogram.  The 
laser  interferogram  is  roughly  a  cosine  wave  which  crosses  zero  at  X  =  j.\/2 

where  X  =  HeNe  and  j  =  1,  2,  3 -  Since  the  laser  operates  in  air,  the 

measured  path  difference  Xm  is  sampled  with  a  spacing  of  AXm  given  by. 


AXm  =  j  •  4AAir(HeNe),  j  =  1,  2...  (7) 

where  AAir  is  the  wavelength  of  the  helium-neon  laser  in  air.  In  the 
calculation  of  Bc(oc),  it  is  assumed  that  the  path  difference  was  sampled 
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From  equation  {12),  equation  (11)  can  be  written: 


Bc(oc>  '  fW  s(STf®iiT  ‘  °c)dV  (13) 

0 

Since  only  the  alteration  of  the  wavenumber  scale  is  of  interest,  inspection 
of  equation  (13)  shows  that  the  delta  function  is  zero  except  at, 

°c  =  n(HeNe)  (14) 

and  since  the  observed  wavenumbers  are  in  air,  aQ  can  be  written, 


a  =  n(o..  ) 

o  v  Vac' 


Vac 


(15) 


where  n(°Vac)  is  the  index  of  refraction  of  air  at  the  vacuum  wavenumber 
of  the  spectrum.  From  equations  (14)  and  (15),  the  calculated  wavenumbers, 
oc,  can  be  related  to  the  vacuum  wavenumber  of  the  observed  spectrum  by, 

(16) 


°c 


n(o»ac> 


n(HeNe)  Vac 


The  ratio  of  the  indicies  of  refraction  of  air  in  equation  (16)  is  less  than 
one  for  aVac  lying  in  the  IR,  so  the  calculated  wavenumbers  are  less  than  the 
true  (vacuum)  wavenumbers.  Consider  the  error.  Act,  resulting  from  equation  (16) 

A a  =  a,,  -  a  =  a,. 

Vac  c  Vac 

At  room  temperature  (22°C)  and  standard  pressure  (760mm)  the  error  in  the 
calculated  wavenumbers  at  2000  cm-1  is  .008  cm-1  and  at  4000  cm"’  the  er>'or 
is  .015  cm'1.  It  is  therefore  necessary  that  this  error  be  considered  if  high 
precision  wavenumber  determination  is  desired.  It  should  be  pointed  out  that 
the  errors  described  by  equations  (1)  and  (16)  lead  to  a  lower  wavenumber 
scale  than  the  true  scale,  thus  their  effects  are  additive. 

A  number  of  other  errors  can  effect  the  precision  of  the  wavenumber 
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scale  on  the  FTS  system;  however,  the  errors  summarized  by  equations  (1)  and 
(16)  are  both  significant  and  readily  correctable.  Problems  such  as  sampling 
errors,  laser  instability,  and  mirror- velocity  fluctuations  can  lead  to  errors 
in  wavenumber  calibration,  but  these  are  normally  characteristics  of  a  particular 
FTS  system,  and  may  not  be  easily  corrected.  While  equation  (16)  can  supply  a 
correction  to  the  calculated  wavenumber  scale  by  knowing  the  index  of  refraction 
of  air,  the  solid  angle  in  equation  (1)  can  be  difficult  quantity  to  accurately 
measure.  One  approach  to  determining  the  solid  angle  term  in  equation  (1)  is 
to  compare  measured  line  positions  with  known  values.  Then,  theoretically,  the 
correction  term  for  equation  (1)  is  known  for  the  remainder  of  the  wave  number 
scale.  This  was  the  scheme  chosen  to  calibrate  the  FTS  wavenumber  scale  resulting 
in  a  calibration  program,  FTCAL,  which  is  described  in  appendix  1 1 1- 1 .  Using 
known  line  positions  (up  to  100),  the  solid  angle  term,  q/4tt  ,  in  equation  (1) 
is  determined  by  the  method  of  least-squares  after  altering  the  known  line 
positions  according  to  equation  (16).  The  reason  for  altering  the  known  line 
positions  before  performing  the  least-squares  calculation  is  so  the  only  error 
remaining  is  that  due  to  equation  (1).  After  the  least-squares  calculation  is 
performed,  corrected  values  for  the  measured  line  positions  are  computed  as  follows 

a  (measured) 

'  1  -  SV4tt 

o  (corrected)  =  aj  •  -  (171 

n\0i  l 

The  corrected  values  are  compared  with  the  (unaltered)  known  positions  to 
determine  the  errors  between  the  two  sets.  If  the  user  has  specified  an 
inclusion/rejection  criterion  (the  parameter  PR  in  the  program),  the  errors 
are  compared  with  the  value  of  PR  to  determine  if  a  line  should  be  excluded 
from  the  least-squares  calculation. 
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If  the  largest  error  exceeds  PR,  the  measured  line  position  which  yielded  that 
error  is  excluded  from  the  set  of  measured  values,  and  the  least-squares  cal¬ 
culation  is  repeated.  This  process  is  repeated  until  only  the  lines  yielding 
an  error  less  than  PR  are  included  in  the  least-squares  calculation  of  Q./ 4ir. 

If  PR  is  set  equal  to  zero,  all  lines  are  included.  After  a  value  for  fi/4-rr 
is  found,  unknown  lines  are  read- in  and  corrected  according  to  equations  (17). 

The  index  of  refraction  of  air  is  calculated  using  Edlen's  formula  and  corrected 
for  non-standard  temperature  and  pressure  using  ideal  gas  law  corrections  (21. 

The  user  of  FTCAL  must  specify  the  temperature  (deg.  C)  and  barometric  pressure 
(nri  of  Hg)  at  which  the  standard  lines  were  measured  and  the  unknown  measurements 
made.  It  should  be  pointed  out  that  if  the  temperature  differs  by  less  than 
±5°C  and  the  pressure  differs  by  less  than  +20  torr  between  the  calibration 
measurement  and  the  unknown  measurement,  the  error  resulting  from  assuming 
both  were  taken  under  the  same  conditions  is  less  than  .001  cm'1  over  the 
wavenumber  range  400  cm'1  to  4000  cm"1.  For  a  detailed  description  of  FTCAL, 
the  reader  is  referred  to  Appendix  I. 

Another  area  that  needed  investigating  before  making  precision  line 
position  measurements  involved  the  best  methods  for  determining  the  line 
center  with  the  FTS  system.  Two  basic  options  are  available  to  the  FTS  user. 
Firstly,  the  user  may  use  the  intensified  cursor  to  directly  read-out  the  line 
positions  from  the  c.r.t.  display  and  secondly,  the  user  may  read  the  line 
position  values  from  a  hard-copy  plot  of  the  data.  The  former  method  is  very 
convenient,  but  suffers  from  poor  consistency.  The  latter  method  is  cum¬ 
bersome  and  time  consuming  but  yields  good  measurement  consistency.  Using 
the  carbon  monoxide  fundamental  as  a  known  set  of  lines  and  the  overtone  as 
an  unknown  set  (both  are  known  to  better  than  .0002  cm"1)  evaluation  of  the 
FTS  system  and  the  FTCAL  program  was  begun.  From  these  measurements,  a 
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number  of  guidelines  developed.  These  are  summarized  below: 

1.  Dispersive  elements  (such  as  lenses)  in  the  collimating  optics  lead 
to  non-linearities  in  equation  (1)  and  invalidate  that  correction  if  the 
unknown  and  known  lines  differ  greatly  in  wavenumber. 

2.  Use  of  the  intensified  cursor  to  determine  line  positions  can  at  best 
yield  a  precision  of  *.005  cm'1. 

3.  Use  of  the  plotting  routine  with  the  FIT  parameter  set  to  YS  can  lead 
to  a  precision  of  better  than  *.002  cm-1.  Use  of  the  FIT  routine  is 
imperative  for  good  interpolation  between  data  points. 

4.  Line  centers  should  be  determined  by  using  paper  with  a  grid  or  over¬ 
laying  plots  on  a  grid.  Plots  should  always  be  made  of  absorbance  files 
unless  absorbance  values  are  small  (<0.I).  Plots  should  be  made  using  a 
scale  factor  of  0.1  cm'1  per  inch  or  less. 

It  is  felt  that  a  precision  of  better  than  ±.002  cm'1  is  obtainable 
over  the  wavenumber  -ange  4 00  cm'1  to  4000  cm'1  and  perhaps  ±.001  cm'1  in 
the  wavenumber  range  below  2000  cm'1. 
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SECTION  II 


ROTATIONAL  AND  VIBRATIONAL  POPULATION 
DISTRIBUTION  IN  CARBON  MONOXIDE 

Determination  of  the  rotational  temperature  and  vibrational  populations 
of  carbon  monoxide  in  a  low  pressure  gas  discharge  has  been  accomplished.  A 
computer  code  has  been  developed  which  allows  the  user  to  input  low  resolution 
emission  data  from  either  the  fundamental,  first,  or  second  overtone  of  carbon 
monoxide  and  outputs  to  the  user  the  rotational  temperature  and  the  fraction 
of  molecules  in  each  vibrational  level.  This  investigation  resulted  in  a 
paper  being  presented  at  a  Mini  Symposium  on  Aerospace  Systems  Technology  - 
Present  and  Future,  sponsored  by  the  American  Institute  of  Aeronautics  and 
Astronautics. 

Figure  4  shows  the  experimental  setup  used  for  making  the  measure¬ 
ments.  Radiation  leaving  the  discharge  tube  passed  through  the  adjustable 
aperture,  through  the  fourier  transform  spectrometer  (FTS ) ,  and  was  then 
focused  on  the  PbSe  detector.  The  discharge  tube  was  a  cylinder  roughly  2.5  cm. 
in.  diameter  and  50  cm.  in  length.  The  entire  optical  path  was  purged  with 
dry  air  which  was  also  free  of  carbon  dioxide.  Both  water  vapor  and  carbon 
dioxide  parasitical ly  absorb  the  infrared  radiation  from  the  discharge.  The 
spectral  regions  in  which  these  molecules  absorb  is  illustrated  in  Figure 
3,  and  partially  overlaps  the  emission  from  CO. 

When  absorption  due  to  CO^  and  H^O  is  not  present,  the  spectrum  of  the 

emission  from  a  CO-N^-He  discharge  is  illustrated  in  Figure  5.  The  band 

at  2000  cm'1  is  due  to  many  overlapping  fundamental  transitions  AV=1  while 

the  band  extending  from  2800  to  4400  cm'1  is  due  to  many  overlapping  first 

overtone  (AV=2)  transitions.  As  shown  schematically  in  Figure  6,  each 

of  the  overtone  transitions  has  a  P  and  an  R  branch,  but  each  successive 
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SCHEMATIC  OF  THE  EXPERIMENT 

FIGURE  4 
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transition  is  shifted  to  lower  wavelength  because  of  vibrational  anharmonicity. 
The  observed  spectrum  is  the  overlap  of  many  slightly  shifted  P  and  R  branches. 
It  should  be  noted  in  Figure  5  that  the  R-branch  of  the  V  =  2  to  V  =  0 
transition  is  not  overlapped  by  any  other  bands,  hence  the  theoretical 
description  of  the  intensity  in  this  region  is  greatly  simplified.  Figure 
7  is  detailed  plot  of  the  CO  first  overtone  emission  showing  the  structure 

resulting  from  more  than  30  vibrational  levels  contributing  to  the  spectrum. 

The  approach  taken  to  the  problem  of  determining  the  rotational  temper¬ 
ature  and  vibrational  populations  from  data  typified  in  Figure  7  was 
guided  by  three  primary  considerations.  First,  determine  the  rotational 
temperature  and  vibrational  populations  that  yield  a  "best  fit"  between  a 
measured  and  a  theoretical  first  overtone  spectrum.  Second,  the  computation 
should  be  efficient  and  free  from  "trial  and  error"  fitting  procedures.  Third, 
because  of  the  non-equii ibrium  nature  of  the  problem,  no  assumption  about 
the  behavior  of  the  vibrational  populations  can  be  made.  In  addition,  two 
assumptions  made  in  the  approach  were  that  the  distribution  of  molecules 
among  rotational  lines  could  be  treated  as  delta  functions. 

The  method  of  solution  developed  from  this  approach  is  outlined  below. 

I .  Correct  the  measured  spectrum  for  the  instrument  response  function. 
The  instrument  response  function  (IRF)  is  found  by  dividing  a 
measured  black  body  spectrum  by  the  theoretical  black  body  spectrum. 
Raw  spectral  data  is  subsequently  corrected  by  dividing  by  the  fRF. 
This  procedure  removes  effects  due  to  beamsplitter  variations, 
detector  sensitivity,  amplifier  gain,  etc.  The  spectra  shown  were 
corrected  in  this  manner. 

II .  Calculate  the  rotational  temperature  from  the  R-branch  of  the 
2-0  transition.  Because  the  2-0  R-branch  is  not  overlapped 
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by  other  bands,  the  shap1  only  depends  on  T 
number  density.  These  parameters  are  determined  by  performing 
a  two  parameter,  non-linear  leas*-  squares  fit  hetween  the  measured 
and  theoretical  2-0  R-branch. 

III.  Calculate  the  vibrational  populations  using  uata  from  full  overtone 
spectrum.  Using  TRqT  from  II,  a  theoretical  overtone  spectrum 
is  generated  which  is  a  function  of  the  number  densities  in  each 
vibrational  level.  A  multiple  linear  regression  analysis  is  then 
used  to  determine  the  set  of  vibrational  populations  that  yield  a 
best  fit. 

IV.  Compare  the  resulting  theoretical  spectrum  with  the  measured  spectrum. 
Step  I  is  performed  on  the  FTS  data  system  while  the  remaining 
steps  are  performed  on  a  CDC  6600  by  the  computer  code  CODIAG.  A 
listing  of  this  code  is  Appendix  II.  Data  points  are  taken  from  the 
FTS  data  system  and  manually  punched  on  cards  for  the  6600. 

A  map  of  the  structure  of  CODIAG  is  shown  in  Figure  S  .  The  first 
step  in  the  analysis  is  to  determine  the  rotational  temperature  from  the 
spectral  data  of  the  R-branch  of  the  2-0  transition.  The  intensity  of  the 
CO  emission  in  this  region  can  be  written  as, 

Ko)  =  N2f(a;TR)  (18) 

where  1(a)  is  the  intensity  at  wavenumber  o  ,  N2  is  the  number  of  CO  mole¬ 
cules  in  the  v= 2  level,  and  f(a;TR)  is  a  function  which  is  non-linear  ir. 
the  rotational  temperature,  TR.  The  behavior  of  the  function  f(a;TR)  is 
the  result  of  the  intensity  distribution  among  the  rotational  lines  of  the 
2-0  R-branch  folded  with  the  spectrometer  instrument  function.  The  unknown 
parameters  N2  and  TR  in  equation  (18)  are  determined  from  a  set  of  measured 
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FIGURE  8  PROGRAM  STRUCTURE  OF  CODIAG 
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intensities  (I(o.)}  using  the  non-linear  least-squares  fitting  procedure 
developed  by  Marquardt  f 4]  and  described  in  Bevington  [5],  The  technique 
is  iterative  and  requires  initial  guesses  for  N  and  TR.  This  procedure 
is  applied  through  the  C0D1AG  SUBROUTINF  ROTFM.  Initial  guesses  and  data 
input  is  through  ROTEM  while  the  actual  least-squares  algorithm  is  contained 
in  SUBROUTINE  CURFIT  which  is  modified  version  of  a  program  appearing  in 
Bevington  [5].  The  right  side  of  equation  (18)  is  calculated  in  FUNCTION  FIO 
which  is  called  by  CURFIT  through  the  intermediary  FUNCTION. 

The  intensity  at  any  point  in  the  calculated  spectrum  can  be  written 


=  ^  Vo(W 


(19) 


■  th  . ... 


where  N.  is  the  number  of  molecules  in  the  j  vibrational  level  and  f.(cs.;TD) 
J  J  i  k 

includes  the  behavior  of  the  rotational  intensity  distribution  folded  with 

the  spectrometer  instrument  function.  Since  TR  has  already  been  found,  write 

f.(a.;T  )  as  f.(o.)'  _  f...  Then  equation  (19)  becomes 
J  '  K  J  *  J» 


!(„,)  =  If  ■  Z 


(20) 


Since  the  measured  intensities  might  contain  a  zero  offset,  the  function 
used  to  describe  the  measured  spectrum  is  written 


I . 
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ao  + 


.f . 

J  J 


(21) 


It  is  desired  then  to  find  the  set  { N j }  from  a  measured  set  of  intensities 
(Y.1.  The  method  chosen  for  this  is  to  minimize  the  function  y2  by  adjusting 
the  set  { N j  } . x?  is  given  by 


(22) 


where  A.  is  the  standard  deviation  of  the  ith  measured  intensity.  At  the 
1 

minimum  of  x2  ♦  the  partial  derivatives  of  x?  with  respect  to  aQ  and  each 

parameter  N.  are  zero  or 
J 

3x2  =3x1=  - *  3x1  =  0 

3a  3N2  3N  (23) 

0  n 

Substitute  equation  (21)  into  equation  (22)  and  apply  condition  (23)  yields 

a  set  of  linear  equations  which  can  be  solved  for  a  and  the  set  {N.}.  The 

o  j 

reader  is  refered  to  Bevington  [5]  for  a  detailed  description  of  the  solution. 

Tha  application  of  this  scheme  to  determine  the  set  { N ^ >  is  through  SUBROUTINE 

VIBDIS  (See  Figure  8).  The  actual  multiple  linear  regression  is  performed 

in  SUBROUTINE  REGRES  which  is  a  modified  version  of  a  program  appearing  in 

Bevington  [5].  The  functions  f..  in  equation  (21)  are  calculated  in  FUNCTION 

J  ■ 

FNV.  In  addition  to  handling  the  data  1.0.  from  the  vibrational  population 

determination,  VIBDIS  also  calculates  values  for  N  and  N  from  an  assumed 

o  i 

Treanor  distribution  6  over  the  lower  vibrational  level.  The  Treanor 
distribution  can  be  written  as 

Ny  =  Nq. Exp  (VO  -  Ey/kT)  (24) 

+  h 

where  Ey  is  the  energy  of  the  V1  vibrational  level  compared  to  the  V=0 
energy,  0  is  a  constant  at  a  fixed  gas  temperature,  k  is  the  Boltzmann 
constant,  and  T  is  the  kinetic  (and  rotational)  temperature  cf  the  gas. 

The  value  of  0  can  be  determined  from  equation  (24)  with  the  previously 
determined  values  of  N2  and  N3  and  is  given  by 

0  =  ln  (  {Jj-4  +  (E,  -  E2)  /kT  (25) 

24 


m 


Using  equation  (25)  and  equation  (24),  along  with  the  calculated  values  for 

and  N3  ,  the  Treanor  values  for  Nq  and  N,  can  be  found. 

As  an  initial  test  of  CODIAG,  data  points  were  taken  from  a  synthetic 

spectrum  generated  from  a  known  T  and  a  known  set  (N.)  and  input  into 

R  J 

CODIAG.  This  test  would  indicate  if  the  program  could  determine  T  and  the 
set  (Nj}from  a  limited  number  of  data  points.  The  results  for  a  30  level 
set  f N j }  were  excellent.  A  more  difficult  test  was  to  determine  the  desired 
parameters  from  a  measured  spectrum.  Data  was  obtained  from  a  liquid  nitrogen 
cooled  C0-N..-He  discharge  with  a  1:2:10  mixture  ratio  and  operating  at  a 
pressure  of  about  1  torr.  The  discharge  length  was  about  50  cm  and  the 
current  and  voltage  were  15  mA  at  4.9  Kv.  In  addition,  about  .15?  of  the  total 
pressure  was  oxygen.  The  measured  overtone  spectrum  is  shown  in  Figure 
as  obtained  with  a  Fourier  Transform  Spectrometer  with  a  resolution  of  about 
10cm"1.  Calibration  of  the  system  was  achieved  by  comparing  a  theoretical 
700°K  black  body  function  with  a  measured  spectrum  from  a  standard  source. 
Eighteen  data  points  were  taken  from  the  R-branch  of  the  2-0  transition  for 
the  determination  of  TR  and  an  additional  seventy  four  data  points  were  taken 
throughout  the  spectrum  for  determination  of  the  vibrational  populations. 

The  data  points  chosen  were  primarily  the  "peaks"  and  "valleys"  of  the 
vibrational  structure  in  the  spectrum.  The  rotational  temperature  was  found 
to  be  207°K  with  an  estimated  uncertainty  of  jd0°K.  The  calculated  vibra¬ 
tional  population  distribution  is  illustrated  in  Figure  9,  which  is  a 
plot  of  the  log  of  the  relative  occupation  number  versus  vibrational 

level.  The  relative  occupation  number  is  given  by  N./ZN.. 

J  J 

Also  shown  in  Figure  9  is  the  population  distribution  from  a 
room  temperature  discharge  shown  in  Figure  10.  The  stars  correspond 
to  the  data  in  Figure  7  and  the  rectangles  correspond  to  the  data  in 
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Figure  10.  The  theoretical  spectra  calculated  from  these  vibrational 
populations  and  rotational  temperature  are  shown  in  Figure  11  and 
Figure  12.  In  these  figures,  the  small  diamonds  indicate  the  data 
points  used  in  the  analysis.  Figure  11  and  Figure  12  can  be  compared 
with  Figure  7  and  Figure  10  as  an  indication  of  the  confidence  in  the 
results  shown  in  Figure  9. 

From  the  results  shown  here,  some  conclusions  can  be  drawn  as  to  the 

« 

limitation  of  the  diagnostics  along  with  some  guidelines  for  determining  its 
sensitivity.  With  a  signal-to-noise  ratio  (S/N)  on  the  order  of  100, 
vibrational  population  data  has  been  calculated  over  5  orders  of  magnitude. 

With  a  S/N  on  the  order  of  65  in  the  2-0  R-branch,  the  estimated  standard 
deviation  in  T  is  less  than  5°/  of  T  .  It  appears  that  a  good  estimate 

of  the  relationship  between  S/N  in  the  R-branch  and  the  standard  deviation 
in  Trqt  is  given  by  (AT/T)  •  (S/N)  =  1.3.  The  chief  sources  of  error  in  the 
calculation  of  the  vibrational  populations  are  probably  systematic.  Errors 
in  correcting  for  system  sensitivity  and  errors  in  calculating  the  CO  transition 
probabilities  are  probably  large  compared  to  the  fluctuations  in  the  vibrational 
populations  due  to  a  S/N  greater  than  20.  Measurement  times  can  be  shortened, 
down  to  the  time  required  for  a  single  scan  of  the  FTS,  with  a  corresponding 
decrease  in  S/N.  Single  scan  times  can  be  less  than  0.1  sec.  with  an  effective 
noise  equivalent  power  of  -  7  x  10" 10  watt/cm"1  at  4000  cm"’.  This  effective 
NEP  includes  detector  noise  plus  electronic  noise  and  is  adjusted  for  the  sys¬ 
tem  response  function.  Flence,  if  the  radiant  intensity  of  a  source  is  known 
and  the  solid  angle  subtended  at  the  FTS  detector,  this  effective  NEP  can  be 
used  to  give  an  estimate  of  the  single  scan  S/N  that  can  be  expected.  Inherent 
in  the  FTS  technique  is  the  need  for  sources  to  subtend  a  small  solid  angle 
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at  the  detector;  however,  for  low  resolution  this  ,-egu i rement  is  not  as 
stringent  as  needed  for  high  resolution.  For  the  resolution  needed  in  this 
analysis,  the  source  may  subtend  a  solid  anqle  as  large  as  6  x  10  stor. 
without  degrading  the  technique.  Lastly,  it  should  be  pointed  out  that 
absorption  by  water  vapor  in  the  optical  path  would  completely  prohibit 
obtaining  vibrational  population  data;  however,  T  is  still  measureable 
because  the  2-0  R-branch  lies  outside  the  water  vapor  band. 
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SECTION  III 


ABSORPTION  BY  DISCHARGE  IMPURITIES 

A  number  of  gas  mixes  which  model  the  gas  found  in  the  electron  beam 
closed  cycle  discharge  were  studied  in  an  effort  to  determine  detection 
sensitivities  for  a  number  of  minor  species. 

The  first  trace  species  studied  were  N20,  NO,  and  N0?.  These  were 
chosen  for  their  strong  absorption  bands  in  the  PbSe  region  and  their 
suspected  presence  in  C02  laser  discharges.  Gas  mixes  of  3  parts  N? ,  1 
part  CO2  and  various  concentrations  of  the  nitrogen  oxides  were  made  up  and 
their  absorption  spectra  obtained.  Total  pressure  in  the  White  cell  was 
kept  at  200  torr.  This  was  done  for  the  following  reasons.  First,  below 
one  half  atmosphere  the  line  widths  of  the  trace  species  are  on  the  order 
of  the  best  spectral  resolution  obtainable  with  the  FTS  System  (.06  cm"'). 

Thus,  below  -  380  torr,  the  observed  line  widths  should  be  pressure  independent 
and  interference  between  neighboring  lines  will  be  at  a  minimum.  Second, 
the  ratio  of  the  volumes  of  the  available  sample  cylinders  to  the  White  cell 
volume  is  about  1:5,  so  with  an  EBCC  pressure  of  1  atm.,  the  White  ceil  pressure 
should  be  on  the  order  of  150  torr.  Running  the  gas  mixes  at  200  torr  was  a 
reasonably  good  approximation  to  the  conditions  under  which  samples  were  to 
be  studied. 

With  the  gas  mix  just  discussed,  N20  shows  two  reasonably  strong 
absorption  bands  that  are  free  from  interfering  species.  The  bands  originate 
from  the  00°1  -*•  00°0  and  20°0  -*■  00°0  transitions  and  are  centered  at  2223.8 
cm-1  and  2563.3  cm-1  respectively  [7].  Figure  13  shows  a  portion  of 
the  P-branch  of  the  00°1  -*■  00°0  band  while  Figure  14  shows  both  branches 
of  the  20°0  -►  00 °0  band.  These  spectra  were  obtained  with  100  ppm  N?0 
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in  the  3:1  N2:C02  gas  mix  previously  discussed.  The  P-branch  of  the 
00°1  «•  00°0  transition  lies  at  the  edge  of  strong  absorption  by  a  number 
of  CO?  bands;  thus  there  is  some  interference  with  the  lower  J-value  N?0 
lines.  The  R-branch  of  the  00°1  -*■  00°0  transition  is  completely  obscured 
by  the  C02  absorption.  Figure  15  shows  the  same  spectral  region  as 
Figure  13  only  the  N20  is  missing  from  the  gas  mix  and  the  CO  pressure 
is  somewhat  higher.  It  can  be  seen  that  above  about  2200  cm'1  interference 
from  C02  prohibits  the  used  of  the  N20  lines  for  determining  its  presence. 
The  20°0  00°0  band  lies  in  an  interference-free  region,  but  lacks  the 

strong  absorption  of  the  00°1  -*•  00°0  band.  Table  lisa  summary  of  the 
results  for  the  two  N20  bands  along  with  expected  band  strengths  for  the 
two  bands  studied.  Lines  are  chosen  based  on  the  criteria  of  being 
relatively  interference-free  and  showing  strong  absorption.  The  effective 
absorption  coefficient  was  calculated  from  the  slope  of  the  straight  line 
found  by  a  least  squares  fit  of  concentration  dependence  of  absorption. 
Errors  were  calculated  assuming  a  1 %  error  in  the  transmission  measurements 
which  was  roughly  the  noise  level  of  2000  sec.  run  with  the  PbSe  detector. 
The  transmission  through  the  gas  is  given  by: 

I  =  I0e"ap  (26) 

where  (in  this  case),  a  is  in  units  of  torr"1  and  hence  p  is  the  partial 
pressure  in  torr.  The  estimated  detection  limit  for  each  band  shown  in 
Table  1  is  found  by  assuming  a  value  of  I / I 0  of  .995  can  be  measured  (with 
the  new  InSb  detector),  and  finding  the  pressure  that  yields  this  value. 

The  estimate  is  good  at  a  total  pressure  of  200  torr  and  changes  with  total 
pressure.  Below  about  on  half  atmosphere,  the  detection  limit  is  expected 
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R(17) 


to  increase  linearly  with  decreasing  total  pressure  while  above  one-half 
atmosphere  the  detection  limit  should  be  roughly  constant  due  to  Dressure 
broadening.  For  example,  for  the  00°1  -  00°0  band  the  minimum  detectable 
amount  of  N20  at  100  torr  total  pressure  would  be  rouqhly  twice  the  value 
given  .  hile  at  atmospheric  pressure  the  value  would  be  somewhat  less  that 
the  given  value. 

The  fundamental  band  of  nitric  oxide,  entered  near  1876  cm-1,  lies  in 
a  region  of  strong  water  vapor  absorption.  While  the  P-branch  lies  in  the 
stronger  part  of  the  H20  band,  the  R-branch  is  over-lapped  by  the  1 1 1 0( 2 )  -> 
00°0  band  of  C02  [7]  (the  number  in  parenthesis  locates  the  level  in  the 
Fermi  resonating  group).  Due  to  the  high  number  densities  of  C02  in  this 
experiment,  it  was  decided  to  use  lines  of  the  P-branch  that  were  not  too 
strongly  overlapped  by  H20  lines.  Table  2  gives  the  results  of  the 
measurements  on  NO  and  Figure  16  shows  a  typical  absorption  spectrum  of 
the  selected  NO  lines.  The  numbers  arrived  at  in  Table  2  were  found  in 
the  same  manner  as  in  the  N20  data  with  one  exception.  The  estimated 
detection  limit  was  based  on  being  able  to  measure  a  transmittance  of  99% 
as  opposed  to  99.5%  for  the  N20  data.  This  was  done  because  of  the  poorer 
signal-to-noise  level  in  the  1800  cm’1  region  of  the  InSb  detector. 

The  strongest  N02  band  in  the  IR  is  the  v3  band  centered  at  1362.03  cm”1 
[Ref.  9] .  An  average  of  the  bandstrength  measurements  for  this  band  yield  a 
value  f 9]  of  848  cm’ '/cm- atm.  which  is  somewhat  larger  than  one  half  the 
N20  v3  bandstrength.  The  difficulty  in  using  this  band  for  detecting  N02 
is  that  it  lies  in  the  middle  of  a  very  strong  water  vapor  band  which 
blocks  most  of  the  N02  signal.  Figure  17  shows  the  water  vapor  absorption  in 

this  region  as  a  result  of  the  spectrometer  operating  in  laboratory  air.  The 
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spectrum  was  obtained  using  the  IlgCdlo  half  of  the  “sandwich"  detector  and 
the  interferometer  fitted  with  the  Fe2  03-CaF?  beamsplitter.  The  cutoff  of 
signal  near  1100  cm  1  is  due  to  the  CaF2  transmission  while  the  cutoff  near 
1850  cm"1  is  due  to  the  transmission  through  the  InSb  part  of  the  "sandwich". 
From  Figure  17,  it  is  clear  that  it  will  be  difficult  to  observe  N02 
absorption  in  this  region;  however,  there  is  a  small  "window"  in  the 
water  vapor  centered  at  1586  cm-1  which  is  shown  in  Figure  18.  Figure  19 
shows  the  absorption  in  this  region  due  to  50  ppm  N02  in  150  torr  N2  and 
indicates  that  it  should  be  possible  to  use  this  region  to  detect  N02.  An 
alternative  band  which  could  be  used  is  the  +v3  band  of  N02  centered  near 
2906  cm'1.  This  band  is  free  of  interfering  species;  however,  the  band 
strength  is  only  17  cm_1/cm-atm  (average  of  measurements  [10])  and  hence 
will  yield  a  higher  detection  limit. 

Since  the  CO  spectrum  is  well  known,  only  mention  of  the  interference 
free  lines  will  be  given  here.  In  the  fundamental  band,  the  P(3),  P(2), 

P(l),  and  R(0)  lines  were  selected  for  being  free  from  interference  by  either 
C02  or  N20  lines  while  in  the  first  overtone,  all  lines  were  considered 
essentially  interference  free. 
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CALIBRATION  PROCEDURE  FOP.  IMPURITY 


CONCENTRATION 


The  easiest  measured  quantity  that  can  be  related  to  the  partial  pressure 
of  a  particular  species  is  the  transmittance  of  the  gas  at  one  of  the  absorbing 
wavelengths.  If  radiation  from  a  source  passes  through  the  gas,  the  trans¬ 
mittance  of  the  gas  can  be  written 

T ( v )  =  100?  x  (27) 

where  I(v)  is  the  intensity  of  the  radiation  reaching  the  detector  after 
passing  through  the  gas,  and  I0(v)  is  the  intensity  of  the  radiation  reaching 
the  detector  when  no  gas  is  present.  Theoretically,  this  can  be  related  to 
the  partial  pressure  of  the  absorber  by 


T(v)  =  100?  x  F(vi-v')  •  exp  -p-x- k(v' )  dv' 


In  equation  (28),  p  is  the  partial  pressure  of  the  absorber,  x  is  the  path- 
length  over  which  the  absorption  occurs,  k(v' )  is  the  absorption  coefficient 
for  the  particular  species  as  a  function  of  wavenumber  (or  wavelength), 
and  F(v-v')  is  the  spectrometer  bandpass  function,  normalized  according  to 


F(v)dv  =  1. 


Ik 


If  the  spectrometer  bandpass  function  is  much  "narrower"  than  the  width  of 
a  spectral  line,  F(v-v')  in  equation  (28)  goes  over  to  a  delta  function  and 
the  measured  transmittance  is  just 


T(v)  =  100?  x  exp  -p-x-k(v) 
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In  this  case,  if  one  knows  the  absorption  coefficient  and  the  pathlength, 
the  partial  pressure  can  be  calculated  directly. 

Unfortunately,  the  situation  that  exists  when  samples  are  analyzed 
with  the  White  cell-FTS  system  is  not  so  simple.  The  gas  samples  are  normally 
at  pressures  of  100  torr  or  less  which  results  in  spectral  linewidths  on  the 
order  of  0.01  cm-1  (FWHM)  or  less,  while  the  best  instrument  resolution  is  on 
or  order  of  0.06  cm-1.  Under  these  conditions  the  functional  dependence 
of  the  transmittance  on  the  partial  pressure  of  the  species  becomes  more  com¬ 
plex.  The  simplest  approach  to  solving  this  problem  is  to  use  calibration 
curves  to  relate  the  measured  transmittance  to  the  partial  pressure  of  the 
species.  In  this  procedure,  a  series  of  spectra  are  taken  of  gas  mixes  which 
contain  known  amounts  of  the  species  of  interest  and  a  calibration  curve  is 
constructed  for  each  absorption  line  of  each  species.  Then  for  a  particular 
absorption  line,  one  can  write 

\  =  Ak(p)  (31) 

where  the  index  k  is  used  to  indicate  a  particular  line  and  is  the  particular 
function  that  describes  the  transmittance  of  that  line.  This  was  the  procedure 
initially  chosen  for  calibration  of  the  White-cell  FTS  system.  The  difficulty 
with  using  the  aforementioned  calibration  procedure  lies  in  the  dependence 
of  the  function  A^  on  the  other  system  parameters.  For  example,  if  the  resolution 
of  the  spectrometer  is  changed,  a  new  set  of  A^'s  must  be  generated  for  each 
new  resolution  setting.  More  important  however,  A^  is  strongly  dependent  on 
the  background  gas  pressure  and  weakly  dependent  on  the  composition  of  the  back¬ 
ground  gas.  Thus  a  separate  set  of  calibration  curves  would  have  to  be  constructed 
for  each  total  pressure  run  in  the  EBCC  system  or  in  the  TEA  laser  experiment, 
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and  for  each  gas  mix  chosen.  For  an  experiment  such  as  the  EBCC  system  or 
the  TEA  laser,  where  one  would  like  to  vary  these  parameters  at  will,  producing 
calibration  data  becomes  exceedingly  time  consuming.  An  alternate  calibration 
procedure,  though  perhaps  less  accurate,  is  to  work  directly  with  equation 
(28)  to  interpolate  and  extrapolate  a  manageable  amount  of  calibration  data. 
This  has  been  the  procedure  adopted  for  the  White  cell-FTS  system. 

To  see  how  equation  (28)  depends  on  the  background  pressure  and  gas 
mix,  one  needs  to  consider  the  absorption  coefficient  k(v).  For  an  isolated 
spectral  line,  k(v)  can  be  written  [11] 

u 

k(v)  =  ~  V(w,y)  (32) 

where  V(w,y)  is  the  Voigt  function  given  by 


exp(-t2)dt 
yr  +  (w-t) 


(33) 


In  equation  (32),  S  is  the  integrated  line  strength  and  aD  is  the  Doppler 
half-width  of  the  line  given  by 


where  \>0  is  the  wavenumber  of  the  line,  k  is  Boltzmann's  constant,  T  is  the 
gas  temperature  in  Kelvin.  M  is  the  mass  of  the  molecule,  and  c  is  the  speed 
of  light.  In  equations  (32)  and  (33),  w  and  y  are  defined  by 


(36) 


where  v  is  the  wavenumber  at  which  k(v)  is  desired,  v0  and  aD  are  as  before, 
and  a  is  the  Lorentz  half-width  of  the  line.  The  Lorentz  width  is  made  up 

Li 

of  two  components,  the  natural  width  aN  and  the  collision  width  ac.  For 
the  range  of  pressures  studied  here,  <<  ac  so  that  for  purposes  of  this 
study  aL  =  ac.  The  collision  half-width  of  the  line  is  given  by 


a 


C 


y  .p. 
T1 


(37) 


where  y.  is  the  line  broadening  coefficient  for  the  line  under  study,  broad¬ 
ened  by  the  ith  component  in  the  gas  mixture  and  p^^  is  the  partial  pressure 
of  the  ith  component.  From  equation  (37)  it  can  be  seen  that  the  collision 
half-width,  and  hence  the  absorption  coefficient  k(v),  depend  on  the  total 
pressure  of  the  gas  and  the  particular  species  present.  To  use  equation  (28) 
to  calculate  partial  pressures  one  must  therefore  know  the  line  strength  S 
of  the  particular  line,  the  molecular  weight  of  the  absorbing  species,  and 
the  broadening  coefficients  for  each  major  species  present  in  the  qas  mix. 

In  addition,  the  partial  pressures  of  the  major  components  of  the  gas  mix 
must  also  be  known. 

Since  most  of  the  data  mentioned  in  the  previous  paragraph  could  be 
obtained  from  the  open  literature,  one  possible  calibration  scheme  would 
consist  of  gathering  the  existing  data  and  calculating  the  partial  pressures 
of  the  major  species  directly.  Since  it  was  desired  to  develop  a  calibration 
scheme  consistent  with  the  calibration  data  that  had  already  been  obtained, 
a  somewhat  different  scheme  was  chosen.  It  was  found  in  the  calculations 
that  errors  in  the  assumed  widths  of  the  lines  were  of  less  importance 
than  errors  in  the  individual  line  strengths.  In  addition,  measurement 

48 


d 


of  broadening  coefficients  would  require  many  additional  calibration  measurements 
which  would  be  quite  time  consuming.  So  for  these  reasons,  the  calibration 
scheme  chosen  involved  using  broadening  coefficients  from  the  open  literature 
along  with  line  strengths  calculated  from  the  calibration  data.  Before  any  of 
this  could  be  implemented,  a  number  of  computer  codes  had  to  be  developed  in 
order  to  evaluate  equation  (28).  To  calculate  the  Voigt  function  in  equation 
(32)  a  code  due  to  Armstrong  [11]  was  modified  slightly  and  implemented.  This 
function  subprogram  is  called  VOIGT.  To  evaluate  equation  (28)  a  code  called 
F0LD1  was  written  which  uses  Simpson's  one-third  rule  to  evaluate  the  integral. 

The  step  size  in  F0LD1  is  automatically  adjusted  to  handle  variations  in  line 
width  and  the  limits  on  the  integral  are  adjusted  in  accordance  with  the 
instrument  resolution.  The  spectrometer  bandpass  function  is  supplied  by  a 
function  subprogram  called  BPASS1  which  was  written  especially  for  this  FTS 
system.  To  calculate  the  strength  of  a  particular  line  it  was  necessary  to 
"solve"  equation  (28)  for  S  (from  equation  (32))  given  a  particular  line  center 
transmittance,  partial  pressure,  and  line  width.  For  this,  a  code  called 
SEEKER2  was  written  which  solves  equation  (28)  for  S  by  an  iterative  procedure. 
Finally,  after  the  line  strengths  were  estimated  a  code  was  needed  to  solve 
equation  (28)  for  p  given  k(v0)  and  T(v0)  for  an  "unknown"  gas  sample.  A  code 
called  SEEKER1  was  written  for  this  purpose  and  like  SEEKER2  uses  an  iterative 
approach.  Data  arrived  at  for  use  in  measuring  minor  species  in  the  White 

cell  FTS  system  is  presented  in  Tables  3  through  7.  The  lines  in  the  tables 
were  primarily  selected  because  they  are  free  from  interference  by  strong  water 
vapor  lines,  hot  bands  of  C02,  and  lines  of  other  minor  species.  Whenever 
possible,  two  bands  of  the  same  molecule  were  used  (weak  and  strong)  to  increase 
the  dynamic  range  of  the  measurement  and  further  reduce  the  chances  of  interference. 
The  line  strengths  given  in  the  third  column  of  each  table  were  arrived  at  by 
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Calculated  from  data  of  Referenced.  Uncertainty  is  probably  less  than  ±5% 
From  Reference  13 
From  Reference  14 
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See  Note  g  of  Table  III- 


first  taking  the  averages  of  the  strengths  calculated  from  each  calibration 
curve  by  SEEKER2,  then  "smoothing"  the  values  to  fit  theoretical  relative 
intensities  at  296°K.  This  "smoothing"  was  done  by  first  calculating  a  band 
strength  from  each  line  strength  using  theoretical  relative  intensities.  Then 
an  average  band  strength  and  its  standard  deviation  were  calculated  and  a  set 
of  line  strength  derived  from  this  average  band  strength  using  the  theoretical 
relative  intensities.  These  final  line  strengths  were  taken  as  the  "smoothed" 
values.  The  data  in  Tables  3  through  7  is  now  on  computer  cards  and 
used  by  a  driving  program  called  TRACER.  TRACER  uses  SEEKER1  and  specific 
line  parameter  data  to  arrive  at  partial  pressures  of  CO,  N20,  and/or  NO  in 
a  given  sample  of  gas.  A  listing  of  TRACER  is  given  in  Appendix  III. 


SECTION  V 


IMPURITIES  IN  CO,  LASER  DISCHARGES 

Gas  samples  from  the  electron  beam  closed  cycle  laser  and  from  a  TEA  laser 
were  analyzed  to  determine  the  extent  of  dissociation  and  gas  chemical  effects 
in  the  discharge.  Before  discussing  the  results  in  general,  the  results  from 
two  gas  samples  will  be  discussed  in  detail. 

Gas  samples  from  the  EBCC  loop  were  always  taken  before  and  after  running 
the  E-beam  to  insure  that  the  gas  was  initially  clean.  Following  irradiation, 
the  samples  were  introduced  into  the  absorption  cell  and  their  spectra  taken. 
Changes  in  the  gas  were  then  noted  by  comparing  the  soectra  of  the  pre-irradiated 
gas  with  the  irradiated  sample.  The  first  irradiated  sample  of  gas  was  from 
approximately  a  two  hour  E-beam  run,  while  the  second  irradiated  sample  was 
taken  after  about  a  seven  hour  run.  In  both  cases,  the  pre-irradiated  samples 
showed  only  the  spectra  expected  from  the  1:2:3  mix.  For  the  two  hour  irradiated 
sample,  the  only  difference  from  the  pre-i rradiated  gas  was  the  Dresence  of 
CO  lines  in  the  irradiated  sample  spectrum.  The  seven  hour  irradiated  sample, 
however,  showed  significant  changes.  For  example,  Figure  20  shows  a  portion 
of  the  absorption  spectrum  of  the  pre-irradiated  gas  while  Figure  21  shows 
the  same  region  of  the  spectrum  of  the  irradiated  sample.  In  the  irradiated 
sample  spectrum  both  N,0  and  CO  lines  can  be  clearly  seen  in  addition  to  the 
CO,  lines  from  the  1:2:3  mix.  Figure  22  shows  a  calibration  spectrum  of 
38  ppm  N,0  in  the  1:2:3  mix  and  indicates  the  order  of  magnitude  of  the  N20 
absorption.  In  addition  to  the  N20  lines  and  the  CO  lines,  lines  due  to  NO 
were  also  seen  in  the  irradiated  sample  which  were  not  seen  in  the  pre- 
irradiated  gas.  Figure  23  shows  a  portion  of  the  irradiated  sample  spectrum 
showing  the  NO  lines,  while  Figure  24  shows  a  typical  calibration  spectrum. 
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FIGURE  21  Absorption  Spectrum 
Irradiated  E-beam 


00*201  00*98 


00*89  00*19  00*f»£ 

33NtUlIW9Ntiyi 


00*Z.I  00* 


WAVENUMBERS 


E 

o 


□□•201  00*58 


00*89  00*15  00  *  *i£ 

OONdlllWSNdyi 


00*Z.I  00* 


oo 


59 


u 

w 

0) 

r— • 

4->  •• 
to  V» 

LUCD 

•r-  CD 

-o  c 

-Jh- 

<Q  •*- 
S-  -J 

s- 

OLr-i 

ZZ 

.  o 

2 

a 

«♦- 

O  C7> 

tnii_ 

C 

E  •*- 

D 

s-  o 
+->  jc~« 

CLUJ 

o  on  i 
<D  E 

UJC-J 

0-0)0 
on  •— 

CD  2 

Q.^1 
c  E 

1  UJ 

O  <0  00 

on  *-« 

lucn  • 
iu  in 

Q-  E  l 

O  fO 

atnuj 

O  0)00 
t/>  -O  ^ 

UJQ-Z 

_o  i  oo 

<  UJr-< 

i—  *-• 
au_i 
r-Z 

CO 

C\J 

□  r-UJ 

1 1 1 

aza 

a: 

cdo*-* 

HD 

CD 

CD  XX 

►— « 

U- 

—cod 

jr 

(D 

fCD 


CM 

-ID 

CO 


□ 

_CO 

03 


CD 

-in 

CD 


CD 

-in 

CD 


.r 

-in 

CD 


an 

CD 


o 

on 

CD 


- 1 — 

00*  *6 


- 1 - 1 - 1 — 

00' 9B  00' 28  00' 9 L 

3aNtJliIWSNtiyi 


- 1 — 

00 '0  L 


CD 

CD 


oo'  hsr 


60 


WAVENUMBERS 


o 


UJUUJ 
Dr^Q 
C2CL'-' 
f-X 
Z^D 
DZr- 
r-1  a 
t-z 
a_a_z 
ccclo  * 
o  mx 
cnocc»- 

QQ«-*£XZ 
CL  U 
X  </) 

ooza 

UJCC*-»tJ 

oca. 

CLQ.UJGC 
oca  axi 
u_  •—•in 
zoxa 

r-*  | — D _ i 


00*001  00*  hS 


00  *  88  00*29  00*  9Z.  00*0Z. 

30Ntill I WGNtiy 1 


03 

oa*^ 


61 


WRVENUMBERS 


Table  8  summarizes  the  results  of  the  two  data  sets  and  includes  estimates 
of  the  concentrations  of  the  species  present.  No  unusual  spectral  features 
were  seen  in  either  of  the  samples. 

Measurements  of  the  selected  trace  constituents  in  TEA  laser  gases  have 
been  obtained.  The  laser  used  is  a  Molectron  model  T250  run  at  about  one-half 
atmosphere  (actually  a  C02  TE  laser).  In  the  experiment,  the  laser  is  run 
sealed  off  for  a  specific  number  of  pulses,  then  a  gas  sample  is  taken  and 
analyzed.  The  limit  on  the  number  of  pulses  that  can  be  run  is  determined  by 
the  onset  of  repeated  arcing  of  the  laser.  Figure  25  shows  the  results 
obtained  for  one  series  of  data  runs  of  the  laser.  The  gas  mix  used  was  made 
up  of  260  torr  of  a  1:2:3  mix  (C02:N2:He)  and  40  torr  of  H2.  Repetition  rate 
for  the  laser  was  kept  low  to  reduce  gas  heating  and  was  about  0.75  Hz. 

Pulse  width  was  roughly  1  usee  and  voltage  was  about  25Kv.  In  all  of  the  TE 
laser  results,  the  predominant  NO*  compound  is  NO,  while  in  most  of  the  EBCC 
data  (at  atmospheric  pressure)  the  predominant  species  is  N20.  The  only 
exception  to  this  was  some  EBCC  data  taken  after  only  one-half  hour  irradiation 
of  the  gas.  In  this  very  short  irradiation,  no  N20  lines  were  seen  while  the 
weaker  NO  lines  are  measureable.  Thus,  it  appears  that  except  for  very  short 
runs  of  the  E-beam,  the  predominant  nitrogen  oxide  is  N20  while  in  all  the 
TE  laser  samples  taken  so  far,  NO  is  predominant.  In  all  cases,  no  N02  could 
be  detected. 

A  systematic  study  was  performed  of  gas  decomposition  in  the  EBCC  discharge 
and  the  effect  of  hydrogen  and  oxygen  as  gas  additives.  The  results  were  pre¬ 
sented  at  the  1979  Gaseous  Electronics  Conference  [24].  For  all  the  results 
presented  here,  the  laser  was  operated  with  250  Torr  of  C02  and  500  Torr  of 
N2  for  a  period  of  one  hour  before  the  gas  was  analyzed. 


ESTIMATES  OF  MINOR  SPECIES  CONCENTRATIONS  IN  EBCC  GAS  SAMPLES 


NUMBER  OF  PUL5E5 


FIGURE  25  MINOR  5PECIE5  BUILD  UP  IN  B  C02  TE 


With  a  pure  N2/C02  gas  mix  the  oxides  of  nitrogen  built  up  to  the  10  ppm 
(parts  per  million)  range  while  the  CO  concentration  was  of  the  order  of  a  few 
hundred  ppm.  On  line  absorption  measurements  with  the  tunable  diode  laser  showed 
that  the  build  up  of  CO  was  very  rapid  when  the  laser  was  initially  turned  on, 
and  reached  a  stable  value  after  a  few  hundred  seconds. 

This  laser  system  has  very  good  vacuum  capabilities,  whereas  most  other 
laser  systems  have  larger  vacuum  leaks.  A  small  amount  of  oxygen  was  added  to 
the  gas  mix  to  simulate  the  more  usual  conditions.  The  results,  shown  in 
Figure  26,  indicate  an  increase  in  both  CO  and  the  oxides  of  nitrogen. 

At  a  pressure  of  five  Torr  added  oxygen,  the  concentration  of  CO,  N20,  N02  and 
NO  are  increased  by  factors  of  8,20,  12,  and  6  respectively  over  the  results 
obtained  without  added  oxygen.  Prior  to  the  experiment,  it  was  expected  that 
the  additional  oxygen  would  react  with  CO  to  reform  C02  lowering  the  CO 
concentration.  Instead  the  oxygen  appears  to  enhance  the  formation  of  CO. 

The  effect  of  adding  hydrogen  was  investigated  since  many  C02  systems 
now  use  gas  mixtures  containing  hydrogen  [25].  As  shown  in  Figure  27, 

NO  and  N02  concentrations  decrease  with  added  hydrogen.  The  behavior  of  N20 
might  be  explained  by  dissociation  of  HN02  or  HN03  formed  in  the  discharge 
from  NO  and  N02.  No  HN03  or  HN02  was  observed  in  the  gas  mix,  but  wouldn't 
be  expected  in  large  quantities  if  their  destruction  to  N20  was  rapid  compared 
to  their  formation.  The  sensitivity  to  HN03  detection  was  estimated  at  10  ppm. 

As  hydrogen  is  added  to  the  gas  mix  the  concentrations  of  both  CO  and 
H?0  increase.  This  may  be  because  the  hydrogen  reacts  with  the  dissociatively 
produced  oxygen  to  prevent  the  back  reaction  of  CO  to  reform  C02. 

Finally  hydrogen  and  oxygen  were  added  together  to  the  laser  gas  mix, 
with  the  results  as  shown  in  Figure  28.  As  before  the  concentration  of  NO 
and  N02  after  one  hour  irradiation  was  smaller  when  hydrogen  was  added  than 
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FIGURE  27 


EFFECT  OF  HYDRQEEN  DN  E.B.C.C.  SYSTEM  ER5 
COMPOSITION  AFTER  i  HOUR  I RRRD  I RT I  DN  '• 
OXYEEN  FREE  INITIRL  ER5  MIX. 


RELATIVE  NUMBER  DENSITY  (PPM) 


FIGURE  28 


EFFECT  OF  HYDROGEN  ON  E.B.C.C.  SYSTEM  GRS 
COMPOSITION  RFTER  I  HOUR  I RRRD I RT I  ON : 

5  TORR  OXYGEN  IN  INITIRL  GR5  MIX. 


r\j 


PRRTIRL  PRESSURE  OF  H2  (TORR) 


when  hydrogen  was  absent.  However  the  concentration  of  NO,  N02  and  N20  is 
greater  for  a  given  amount  of  hydrogen  when  both  hydrogen  and  oxygen  were 
added  than  when  hydrogen  alone  was  added.  The  CO  concentration  drops  to  a 
level  below  that  observed  in  the  basic  N2:C02  mix.  To  explain  the  behavior 
of  CO  two  competing  effects  must  be  considered.  The  presence  of  H2  lowers  the 
concentration  of  discharge  produced  oxygen  and  tends  to  increase  the  CO  level 
by  preventing  its  back  reaction  of  reform  C02.  It  also  appears  that  H20  tends 
to  reduce  the  CO  concentration.  The  catalytic  reaction  of  H20  with  CO  to  form 
C02  has  previously  been  observed.  As  in  the  previous  case,  CH4  and  HCN  appear, 
but  only  at  higher  H2  concentrations. 

In  conclusion,  the  addition  of  oxygen  to  the  gas  mix  greatly  increases 
the  production  of  nitrogen  oxides  and  carbon  monoxide.  The  addition  of  hydrogen 
to  an  oxygen  free  gas  mix  produces  water  at  ^lOT'.  of  the  H2  concentration, 
increases  the  production  of  N20  and  CO  while  decreasing  NO  and  N02. 

The  addition  of  hydrogen  and  a  fixed  quantity  of  oxygen  to  the  gas  mix 
produces  water  at  55%  of  the  original  H2  concentration  (up  to  twice  the  added 
02  concentration),  increases  the  N20  production  but  decreases  that  of  CO,  NO, 
and  N02.  When  hydrogen  was  added  to  the  mix,  with  or  without  oxygen,  CH4 
and  HCN  were  produced,  but  no  HN03  or  HN02  was  detected  in  the  irradiated 
gas  mix. 
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APPENDIX  I 


Description  of  the  Fourier  Transfrom  Spectrometer  Calibration  program 
FTCAL.  This  program  is  a  modified  version  of  «  program  appearing  in  Tabjes 
of  Wave  numbers  for  the  Calibration  of  Infrared  Spectrometers  compiled  by 
A.R.H.  Cole,  2nd  edition,  Pergamon  Press,  1977,  pp.  202-205. 

The  program  fits  a  set  of  up  to  100  standard  lines  to  one  of  two  fitting 
functions.  In  normal  usage,  the  standard  lines  are  fit  to  the  function, 

a(standard)  =  Oj  •  -  (A-l) 


where. 


o (measured) 

°i  r=a  L 


( A-2 ) 


and  the  parameter  'a'  is  determined  by  the  method  of  least-squares.  n(HeNe) 
and  n(oj)  are  the  indices  of  refraction  of  air  at  the  wavenumber  of  the  helium- 
neon  laser  and  wavenumber  o^.  The  program  may  also  be  used  to  fit  the 
standard  lines  to  a  non-linear  function  of  the  form, 

2  3 

Oi  =  aj  +  a2o  +  a3o  +  auo  (A- 3) 


where  o  is  the  measured  line  position  and  the  a^'s  are  determined  by  least- 
squares.  The  fitting  function  used  is  determined  by  the  user  by  setting 
the  parameter  SWITCH  to  either  1  or  0.  After  the  fitting  functions  is 
determined,  measured  line  positions  of  unknown  lines  are  read-in  and  corrected 
according  to  either  equation  (A-2)  or  (A-3)  and  equation  (A-l).  For  extra¬ 
polation  out  of  the  wavenumber  region  of  the  standard  lines,  equation  (A-3) 
(SWITCH  =  0)  should  not  be  used.  For  interpolation  within  the  wavenumber 
region  on  the  standard  lines,  either  equation  (A-2)  or  (A-3)  may  be  used. 


Description  of  the  input  parameters: 

PR  -  Sets  the  maximum  allowable  error  (in  cm-1)  between  the  corrected 
measurements  and  the  known  values  of  the  standard  lines  used.  After  a  least- 
squares  fit  is  made  to  the  known  standards,  the  maximum  error  is  found  between 
the  corrected  measurements  and  the  known  line  positions  of  the  standards.  If 
this  maximum  error  exceeds  the  value  PR,  the  measurement  that  yielded  this 
error  is  removed  from  the  set  of  measured  values  and  the  least-squares  fit 
is  repeated.  This  rejection  process  continues  until  the  maximum  error  no 
longer  exceeds  PR.  If  PR  =  0,  all  measured  values  are  used. 

SWITCH  -  This  parameter  determines  which  fitting  function  is  chosen. 

If  SWITCH  =  1,  equation  (A-2),  which  is  linear  in  the  measured  values,  is  used. 

If  SWITCH  =  0,  equation  (A-3),  which  is  non-linear  in  the  measured  values, 
is  used. 

TEMPC  and  PRESC  -  The  air  temperature  (deg.  C)  and  barometric  pressure 
(mm  of  Hg)  at  which  the  standard  lines  were  measured. 

TEMPU  and  PRESU  -  The  air  temperature  (deg.  C)  and  barometric  pressure 
(mm  of  Hg)  at  which  the  standard  lines  were  measured. 

Some  additional  notes  to  the  user  of  FTCAL  are  in  order.  All  cards  read 
by  FTCAL  except  the  first  card,  are  read  in  the  list  directed  read  format  of 
the  CDC6600.  Data  may  be  p iaced  anywhere  on  a  card  and  if  more  than  one  para¬ 
meter  is  to  be  read  from  a  card,  the  parameters  may  be  separated  by  spaces, 
commas,  or  other  symbols  (see  CDC6600  manual).  The  first  number  on  a  standard 
line  card  and  an  unknown  line  card  are  for  line  indentification.  For  standard 
lines  this  can  be  an  integer  from  0  to  999;  for  unknown  lines  this  can  be  an 
integer  from  0  to  9999.  If  the  number  of  standard  lines  used  is  less  than 
that  needed  for  a  least-squares  fit  to  the  fitting  function  chosen  (SWITCH  =  0), 
the  function  is  truncated  starting  with  the  highest  powers  of  the  measured  values. 
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A  listing  of  FTCAL  is  given  on  the  following  pages.  The  comment  cards 
explain  the  order  that  the  data  should  appear  on  cards. 
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APPENDIX  II 


On  the  following  pages  is  a  listing  of  the  program  CODIAG  and  its  associ¬ 
ated  subroutines  and  function  subprograms.  The  programs  are  well  coirmented 
so  little  explanation  will  be  given  here.  One  important  point  is  worth 
mentioning.  The  parameter  NDP  (fourth  card,  first  number)  determines  the 
assumed  resolution  of  the  interferometer.  It  should  reflect  any  effects 
due  to  finite  source  size  and  should  only  include  data  points  on  one  side 
of  the  interferogram. 
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